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We report an experimental and theoretical study of the valence electron spectrum of the dilute 
magnetic semiconductor (DMS) material, Mn doped GaAs. The experimental data are obtained 
through the differences between off- and on-resonance photoemission data for several concentra- 
tions of Mn atoms. The theoretical electronic structure is calculated by means of a combination 
of density-functional theory in the local density approximation and dynamical mean-field theory 
(LDA+DMFT), using exact diagonalisation as an impurity solver. Theory is found to accurately 
reproduce measured data, and illustrates the importance of correlation effects in Mn doped GaAs, 
and that quite likely correlations effects are significant in DMS materials in general. In addition, 
our results demonstrate the presence of Mn states at the top of the valence band, and we argue 
that the double-exchange interaction is a significant driving mechanism for the ferromagnetic order 
in this material. 
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The field of diluted magnetic semiconductors (DMS) 
[U [2] has for the past decade broadened to include sev- 
eral aspects of materials science, e.g. materials synthe- 
sis in bulk and thin film form, structural and magnetic 
characterization, photoelectron spectroscopy for charac- 
terization of the electronic structure, and theoretical ap- 
proaches based on model Hamiltonians or first principles 
calculations. Several reviews have lately been authored, 
e.g. in Refs. |3H21 showing the evolution of the field from 
the first reports on magnetic impurities in semiconduc- 
tors [BJ. A milestone event was the report on carrier- 
mediated ferromagnetism of Mn-doped InAs and GaAs 
P] , where the Mn atoms both act as acceptors and mag- 
netic impurities. 

In practical applications of spintronics, an ordering 
temperature above room temperature is needed, and a 
few systems were predicted [2] to be promising from that 
point of view, namely Mn-doped ZnO and Mn-doped 
GaN. Unfortunately, it has proven very difficult to de- 
termine their ordering temperatures experimentally, and 
values ranging from close to 1000K to just above K 
have been reported. This vast spread of values can pos- 
sibly be linked to difficulties in synthesizing samples with 
a homogeneous distribution of substitutional Mn atoms 
in the host lattice. Mn doped GaAs offers a considerably 
more robust situation from a synthesis point of view, and 
substitutional Mn doping at the Ga sites can reach con- 
centrations up to 6-8 %. For the highest concentrations, 
ordering temperatures (Tq) as high as 170 K [7] to 190 



K [5] have been reported and are now generally accepted. 

On the theory side, it has been suggested [2] that the 
ferromagnetic interaction in Mn doped GaAs is described 
by Vonsovsky-Zener's p-d exchange model [9!. The main 
assumption of this model is that the Mn-3d electrons are 
localized and have energies below or close to the bottom 
of the valence band, which also allows them to induce 
a hole in the host material. This results in a valence 
band moment which is coupled anti-fcrromagnetically to 
the impurity 3d moment. The effective exchange inter- 
actions between the localized moments is mediated by 
the conduction carriers of the host, and is argued to be 
long ranged [2]. Therefore even for dilute doping con- 
centrations, with Mn atoms distributed far apart, a sig- 
nificant ferromagnetic coupling between transition metal 
moments was suggested to be possible (2). 

However, recent experimental studies of Mn doped 
GaAs pTO] point to that the Fermi level is situated es- 
sentially in the band gap, not in the valence band, and 
this casts doubts on the applicability of the p-d exchange 
model. Furthermore, the picture of long-ranged interac- 
tions has also been challenged by calculations based on 
first principles theory. In this case a much more short 
ranged interaction is observed in general, with magnetic 
percolation as a key ingredient, and a decent agreement 
is found between calculated and observed ordering tem- 
perature |TT1 [12] . In addition, the works based on first 
principles theory suggest a different location of the Mn-3d 
states, lying at the top of the valence band, and form- 
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ing a strongly hybridized state |3] . This has implications 
for the mechanism of magnetic ordering, and in Refs. 0] 
and H3l it was suggested that the electronic structure of 
Mn doped GaAs has a significant contribution from the 
double exchange mechanism, and that the range of the in- 
teraction was short. Most of these calculations are based 
on a mean-field like description of the electronic structure 
via density functional theory (DFT) in the local density 
approximation (LDA) or the generalized gradient approx- 
imation (GGA). The main criticism to this theory is that 
a mean-field like description may not be sufficiently ac- 
curate for describing the Mn-3d states, and that correla- 
tion effects due to the strong Coulomb interaction must 
be considered [5]. 

The discussion above indicates that a thorough analy- 
sis of the electronic structure of Mn doped GaAs is highly 
needed, combining experimental data with accurate the- 
ory, to ease the interpretation. Recent studies regard- 
ing the coupling of the electronic structure to the na- 
ture of the magnetic interactions further highlight this 
need [H]. The purpose of this letter is indeed to pro- 
vide such an analysis, and we present here the spectral 
properties of Mn doped GaAs as given by differences be- 
tween off- and on-resonance photoemission data for 0.1, 
1 and 6 % Mn doping. The experimental results are 
compared with electronic structure calculations based 
on a combination of density-functional theory in local 
density approximation and dynamical mean-field theory 
(LDA+DMFT) plfTp], 

The photoemission data were obtained at the Swedish 
synchrotron radiation facility MAX-lab using the surface 
end station of the 1511 undulator beamline [17] . The 
samples were prepared in a local molecular beam epi- 
taxy (MBE) system, and were transferred to the pho- 
toemission station in a portable UHV chamber without 
being exposed to atmosphere. The Mn concentrations 
were determined during growth by means of reflection 
high-energy electron diffraction (RHEED) oscillations, as 
described earlier [T5] . All spectra presented here were ob- 
tained from as-grown samples, i.e. samples not subjected 
to post-growth annealing. 

Valence band spectra were recorded at two photon en- 
ergies, one on the 2p 3 / 2 absorption peak (A in the inset 
of Fig. JTh , the other 2 eV below the peak (B in the inset 
of Fig. [ij , and the Mn-3d resonant emission was obtained 
from the difference between the two spectra. In a pre- 
vious study pjj] corresponding structures were extracted 
by combining spectra from (Ga,Mn)As and clean GaAs 
obtained at the same excitation energy. The main advan- 
tage of the procedure adopted here is that it avoids arti- 
facts that inevitably arise due to different surface states 
on different surfaces. In the extraction of difference spec- 
tra it is of course important, particularly in the most di- 
lute cases, that the spectra are properly normalized and 
aligned in energy. This was accomplished by including 
the Ga-3d peak in the recorded energy range. By means 
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FIG. 1. (color online) Resonant photoemission spectra of Mn- 
3d states in Mn doped GaAs at several concentrations. In 
the inset the Mn 2p 3/2 XAS spectra with 0.3 % and 6 % Mn 
(dashed and solid lines, respectively) are shown. "A" and 
"B" indicate the two photon energies used for extracting the 
resonant photoemission spectra. 



of this peak the spectra could be aligned with a precision 
of 1 meV, and the normalization was checked via the ab- 
sence of any systematic structures in the energy region 
of Ga-3d emission. For the most dilute cases it turned 
out that thermal effects, mainly associated with gradual 
monochromator heating had to be taken into account. To 
minimize these effects, the recording time for each pair 
of spectra was reduced to about 2 minutes. To achieve 
reasonable statistics, a number of difference spectra (typ- 
ically 50-100) were added after proper alignment. 

The theoretical spectra have been calculated by means 
of a LDA+DMFT code [H] based on a full-potential 
linear-muffin tin orbital method (FP-LMTO) gT], which 
is particularly suitable for large supercell geometries. 
The Mn doped GaAs was modeled with a supercell 
Ga26MnAs27, corresponding to a Mn concentration of 
3.75%. The experimental lattice constant a — 10.681 a.u. 
was used, and the Brillouin zone was sampled with a 
mesh of 256 k-points. The LDA+DMFT scheme was 
applied to explicitly treat the local Coulomb interac- 
tion between the localized Mn-3d electrons. The 4-index 
rotationally-invariant Coulomb interaction matrix was 
generated from the Slater parameters F°, F 2 and F 4 . The 
choice of the average Coulomb repulsion F° , which cor- 
responds to the Hubbard U, is rather problematic, since 
no calculations based on constrained LDA or RPA meth- 
ods are found in literature. Therefore we have considered 
values between 4 eV and 7 eV, which are the accepted 
strengths of the Coulomb repulsion for bulk metallic 7- 
Mn [22] and MnO [23] . The main results of the paper are 
presented for the intermediate value U = 6 eV, while the 
effects due to the variation of the Hubbard parameter are 
discussed at the end. F 2 and F 4 are easier to evaluate, 
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FIG. 2. (color online) Resonant photoemission spectrum 
of Mn-3d states in Mn doped GaAs at 1% concentration 
compared to the Mn-3d projected density of states from 
LDA+DMFT, with and without Lorentzian smearing. The 
Fermi level is at zero energy, and no artificial shifts have been 
applied between theoretical and experimental curves. In the 
inset the ferromagnetic (FM) and non-magnetic (NM) LSDA 
results are shown, together with the LDA+U results. 

and therefore were calculated directly from the electronic 
density, as for MnO in Ref. [2H These values correspond 
to the average Hund's exchange parameter J ~ 1 eV. 

The effective impurity problem arising in LDA+DMFT 
has been solved through exact diagonalisation (ED) 
method, as described in Ref. [Ml The fermionic bath 
interacting with the atomic impurity has been approxi- 
mated by means of 22 auxiliary bath sites: 18 bath sites 
were coupled to the strongly hybridizing Mn t 2g states, 
while only 4 bath sites were coupled to the weakly hy- 
bridizing Mn e g states. Notice that the full diagonal- 
isation of a model containing a total of 32 states (22 
bath sites plus 10 atomic orbitals) is already extremely 
demanding from the computational point of view, due 
to the fact that full temperature dependence is consid- 
ered [23]. Finally we should mention that all the calcu- 
lations have been made for the paramagnetic phase at 
T = 400 K and using 1200 fermionic Matsubara frequen- 
cies. The double counting problem has been considered 
in the fully localized limit (FLL) [T5] . 

In Fig. [I] the experimental data for the Mn-3d one- 
particle excitation spectrum are reported for samples 
with a Mn concentration of 0.1%, 1%, and 6%. The 
curves were extracted by taking differences between off- 
and on-resonance photon energies for the Mn-3d cross 
section. In Fig.[2]the experimental spectrum for the sam- 
ple with 1% Mn concentration is directly compared with 
the LDA+DMFT projected density of states (PDOS) of 
the Mn-3d orbitals, for 3.7% Mn concentration. For such 
concentrations of Mn, the direct wave-function overlap 



between different Mn atoms is vanishingly small, and it 
is expected that the spectral properties are rather insen- 
sitive to smaller modifications of the Mn concentration 
(as Fig. [I] also shows) . Hence it is possible to compare 
experiment for a Mn concentration of 1% with theory for 
a 3.7% Mn concentration. Note that we in Fig. [2] show 
theoretical data, unbroadened and directly as obtained 
from the DMFT calculation, as well as theory which has 
been broadened due to life-time effects. The broaden- 
ing considered was obtained from Fermi-Liquid theory, in 
which the life-time broadening grows quadratically with 
the binding energy. As is clear from Fig. [2j both experi- 
ment and theory position the Mn-3d states primarily in 
the valence band of the GaAs host, with an important 
feature that the Mn-3d states extend over a broad range 
of energies. The Fermi level is pinned by a Mn-3d band, 
which is situated at the top of the valence band of the 
GaAs host. This feature is in good agreement with the 
recent experimental results of Ref. I10| and has important 
consequences on the nature of the magnetic interactions, 
as discussed below. The experimental spectrum has a 
main peak at ~ 3.2 eV biding energy, an observation that 
is captured by theory. Furthermore, the observed broad 
shoulder, which is located between 4 and 8 eV binding 
energy, is also well described in the LDA+DMFT exci- 
tation spectrum. Overall, the agreement between theory 
and experiment must be judged as being good. Only in 
the closest region to the Fermi level can one notice some 
discrepancies between theory and experiment, where the 
small peak at ~ 1 eV binding energy has a slightly too 
large intensity in the calculations, when compared to ob- 
servations. 

In the inset of Fig. [2] the Mn-3d PDOS for non- 
magnetic and magnetic LSDA are reported, and we note 
that the agreement with the experimental data is less sat- 
isfactory. The entire spectral weight is in this case con- 
centrated relatively close to the Fermi level, and while 
the position of the main peak is almost correct, the large 
shoulder at higher binding energy is totally absent. If 
static correlations are added [55] , as for the LDA+U sim- 
ulations shown in the same inset, the agreement with 
the experiment does not improve much: from the basi- 
cally correct LDA position the main peak is moved to 
lower energies, depending on the value of U, but it is 
not possible to describe the broad distribution of spec- 
tral weight from the Fermi level down to -8 eV, and in 
particular the observed spectral features in at binding 
energies between 4 and 8 eV. The LDA and LDA+U de- 
ficiencies are related to the fact that they are basically 
single Slater determinant theories, and only the inclusion 
of dynamical many-body effects is able to capture the 
satellite-like features of the Mn-3d excitation spectrum. 
These problematic aspects of correlated electronic struc- 
tures are closely related to recent findings [24 , for the 
charge-transfer insulator NiO, where experimental data 
show a broad distribution of Ni-3d spectral weight not 
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FIG. 3. (color online) LDA+DMFT PDOS of the Mn-3d elec- 
trons and As-4p electrons, close and far from the Mn ions. 
The Fermi level is at zero energy. 



describable within one-particle approximations, but well 
reproduced by LDA+DMFT. 

In Fig. [3]we report the calculated PDOS for the As-4p 
states, which are the main contributors to the valence 
band, both for an As atom far away from the impurity 
Mn atom and for an As atom close to the Mn atom. As 
seen by inspection of the states projected on the As atom 
far from the Mn impurity, the Fermi level is located above 
the band edge of the valence band. The p-states of the As 
atom close to the Mn atom show some hybridization with 
the Mn-3d states. These results rule out the presence of 
a hole in the valence band that extends over a long range, 
it is in fact primarily located on As atoms in the nearest 
vicinity of the Mn atom. 

Our results draw a physical picture of an impurity 
band formation within the band gap, composed of As 
4p states and Mn 3d states, as also suggested by re- 
cent experimental results [TU]. Previous DMFT studies 
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FIG. 4. (color online) LDA+DMFT PDOS of the Mn-3d elec- 
trons for different values of the exchange parameter J. In the 
inset the same information is shown for different values of U. 
The Fermi level is at zero energy. 



on model systems [37] are also consistent with this sce- 
nario, if the parametrized exchange interactions between 
Mn spins and carriers spins are strong enough. Over- 
all our theoretical and experimental electronic structure 
rules out a pure pd-exchange mechanism in Mn doped 
GaAs, and points to the fact that a significant contri- 
bution from the double-exchange mechanism must take 
place [12], where conduction and magnetic electrons be- 
long to the same band. In this situation (described by a 
Hubbard-like model) the exchange interactions can be of 
essentially non-Heisenberg character [THl [2H] . 

We also want to stress that the precise choice of 
the Hubbard parameter U within a reasonable range of 
values, does not lead to qualitative differences in the 
LDA+DMFT spectral properties. In the inset of Fig. [3] 
the Mn-3d PDOS for U=5 eV and U=7 eV are reported. 
Even for these important changes of U, the main features 
of the two curves remain similar to what observed in 
Fig. [2] and only a small redistribution of spectral weight 
at higher binding energies can be observed. In Fig. [4] 
the LDA+DMFT Mn-3d PDOS for U=6 eV and J=0.8 
eV [2S] is compared with our calculated data for J=l eV. 
The general features of the spectrum are the same, but 
the two main peaks are shifted towards the Fermi level in 
the simulation with lower J. From these results we can see 
that a change of the parameters in the calculation only 
worsens the agreement with experimental data and does 
not change the physical picture outlined in this letter. 

In conclusions we have studied the spectral properties 
of Mn doped GaAs, in a joint experimental and theo- 
retical effort, based on LDA+DMFT, using exact diag- 
onalisation as an impurity solver. Theory is found to 
accurately reproduce measured data, which points to the 
importance of correlation effects in Mn doped GaAs, and 
that quite likely correlation effects are significant in DMS 
materials in general. In addition, our results demonstrate 
the presence of Mn states at the top of the valence band, 
which is consistent with recent experimental observations 
[10) . This rules out a pure pd-exchange mechanism, and 
points to that the double-exchange interaction is a signif- 
icant driving mechanism for the ferromagnetic order in 
this material. 
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